Preparing projected entangled pair states on a quantum computer 
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We present a quantum algorithm to prepare injective PEPS on a quantum computer, a class of 
open tensor networks representing quantum states. The run-time of our algorithm scales polyno- 
mially with the inverse of the minimum condition number of the PEPS projectors and, essentially, 
with the inverse of the spectral gap of the PEPS' parent Hamiltonian. 



Projected Entangled Pair States, or PEPS [2], have 
been proposed as a class of quantum states especially 
suited to describe the ground states of local Hamiltoni- 
ans in quantum many-body physics. PEPS are a higher- 
dimensional generalization of the one-dimensional Matrix 
Product States [3], or MPS, for which many interesting 
properties have been proven: For example, MPS prov- 
ably approximate the ground state of ID local Hamilto- 
nians with constant spectral gap [4, 5], exhibit an area 
law [5] as well as an exponential decay of two-point cor- 
relation functions. Furthermore, for each MPS with the 
injectivity property [6] , a parent Hamiltonian can be con- 
structed with this MPS as its unique ground state. MPS 
can also be prepared efficiently on a quantum computer 
[ ]. PEPS however form a much richer class of states, 
and can e.g. represent critical systems and systems with 
topological quantum order [ ]. It is conjectured that all 
ground states of gapped local Hamiltonians in higher di- 
mensions can be represented faithfully as PEPS, and al- 
though there are strong indications for this fact, this has 
not been proven. What is clear, however, is the fact that 
one can also construct parent Hamiltonians for them [6], 
and the PEPS will be the unique ground states of those 
Hamiltonians if the PEPS obeys the so-called injectivity 
condition [ ] . Many physically relevant classes of PEPS 
on lattices are known to be almost always injective, in- 
cluding e.g. the 2D AKLT state [6]. A particularly in- 
teresting subclass of PEPS is the one that consists of all 
those states whose parent Hamiltonian have a gap that 
scales at most as an inverse polynomial as a function of 
the system size: in that local observable (i.e. the 

local Hamiltonian) allows to distinguish the state from all 
other ones, as the ground state always has energy zero 
by construction. It was an open problem [] ] whether 
such states could however be even created on a quantum 
computer, as an algorithm that would allow to prepare 
any PEPS would allow for the solution of PP-complete 
problems [ ]. 

In this article we show how well-conditioned injective 
PEPS can be prepared on a quantum computer effi- 
ciently. The key idea of our approach is to grow the PEPS 
step by step. We demand that not only our final PEPS is 
the unique ground state of its parent local Hamiltonian, 
but also that there exists a sequence of partial sums of 
the local terms of the parent Hamiltonian, such that each 
partial sum has a unique ground state of its own. Based 
on this assumption, the algorithm starts with a physical 



realization of the valence bond pairs as its initial state 
and iteratively performs entangling measurements on the 
virtual particles to map virtual degrees of freedom to 
physical ones, just as in the definition of the PEPS. The 
PEPS is called injective, iff this map is (left) invertible 
which can only be the case if the dimension of the phys- 
ical space is actually at least as large as the dimension 
of the virtual space at each vertex. Preparing a PEPS 
by measurements may seem to require post-selection to 
project onto the right measurement outcome. To over- 
come this issue we use the Marriott- Watrous trick [9, 10] 
of undoing a measurement based on Jordan's lemma [11] 
and combine it with the uniqueness property of injective 
PEPS [6] to prepare the required eigenstates. A key el- 
ement that contributes to the success of this algorithm 
is the fact that the measurements are not done locally, 
such as in the framework of dissipative quantum state 
engineering [12], but globally by running a phase esti- 
mation algorithm that singles out the ground subspace; 
a similar approach was used in the context of the quan- 
tum Metropolis sampling algorithm [ ]. Alternatively, 
methods for eigenpath traversal [17, 18] can also be ap- 
plied [21]. 

Definitions and Results. — Before stating the result, we 
review the definition of PEPS and their essential proper- 
ties. Recall [2, (>] that PEPS are quantum states de- 
fined over an arbitrary graph G = (V, E) such that 
quantum systems of local dimension d are assigned to 
each vertex. We construct the PEPS by assigning to 
each edge e <E E a maximally entangled state Y^iLi N)- 
In this way, a vertex v £ V with degree k gets asso- 
ciated with k virtual D-dimensional systems. Finally, 
a map : € D <g) € D <g) • • • € D i-> € d is applied 

to each vertex, taking the k virtual D-dimensional sys- 
tems to a single physical d-dimensional system. The 
linear map AM is usually called the PEPS projec- 
tor" and is parameterized by tensors A^ as follows: 
= ZUEl... ljk= i^]\i){j 1 ,...,j k \ where A^ is 
a tensor with k indices. The PEPS can now be written 
as \ip) = J2i 1 ,...,t„=i C [{ A i V J}v}\ii, .-.,in} where C means 

the contraction of all tensors A^ according to the edges 
of the graph. In the most general case the virtual index 
dimension D as well as the physical index dimension d 
may also depend on the edges e and vertices v of the in- 
teraction graph, but we suppress this detail in favor of 
simplicity. Note, that w.l.o.g. A<~ v) > may be assumed, 
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since for arbitrary we can choose a local basis by 
performing a polar decomposition, i.e. — U^A^ 

with unitary and A™ > 0. 

A PEPS is called injective [6], if each PEPS pro- 
jector A^ has a left inverse. For some PEPS this may 
only be true, after some local contractions of a constant 
number of PEPS tensors A^ according to the interac- 
tion graph of the PEPS forming new projectors A^ for 
which the condition above holds. Since this blocking can 
be performed efficiently for constant degree graphs, we 
may assume for the remainder of this paper, that it has 
already been performed, such that each individual A^ 
in our input is already injective by itself. Note, that the 
existence of a left inverse allows us to strengthen the as- 
sumption AW > w.l.o.g. to A^ > for all v. 

For injective PEPS, there is a simple construction [6] 
of a 2-local parent Hamiltonian, such that the injective 
PEPS is its unique, zero-energy ground state. This con- 
struction gives a parent Hamiltonian for a quantum sys- 
tem consisting of n particles with d-dimcnsional Hilbcrt 
spaces. 

Let H be a Hermitian matrix with Ao < Ai its smallest 
and second smallest eigenvalues. Then we call A (if) = 
Ai — Ao the spectral gap of H . For any matrix A, the con- 
dition number k(A) is defined as n(A) = ^ max [^j , where 
Cmax(A) and a m i n (A) are the largest and smallest singu- 
lar values of A, respectively. We are now in a position 
to state the performance of our algorithm as our main 
theorem: 

Theorem 1. Let G = (V, E) be an interaction graph 
with bounded degree and some total order defined on V . 
Let {A^} v< zy m be a set of injective PEPS projectors of 
dimension d x D k associated with each vinVup to ver- 
tex t (according to the total vertex order) describing a 
sequence of PEPS \ipt), and let k = max n(A^) be the 

v£LV 

largest condition number of all PEPS projectors. Let A = 
min t A(H t ), where A(H t ) is the spectral gap of the parent 
Hamiltonian H t of the PEPS \ipt) ■ Then there exists a 
quantum algorithm generating the final PEPS \ ip\v\) with 

probability at least 1 — e in time 0{- — h |V^fcer). 

Algorithm. — Conceptually, PEPS are constructed by 
first preparing entangled pair states = X^j N) f° r 
each edge of the interaction graph describing the PEPS, 
and then projecting the k virtual indices associated with 
each vertex to a single physical index. While this con- 
struction is usually considered only a theoretical device, 
the proposed algorithm is indeed simulating the above 
construction for the case of injective PEPS with gapped 
Hamiltonians. This entails making the virtual indices 
physical as well. 

Figure 1 presents our algorithm in pseudo-code. We 
proceed by explaining each step in detail. PEPS con- 
struction starts in step 2 by distributing maximally en- 
tangled states of the desired bond dimension according 
to the interaction graph G — (E,V). The resulting sys- 
tem is the zero-energy ground state of a simple Hamil- 



Input: Interaction graph G — (V,E) with degree bound k 
and total vertex order. For each v £ V , d x D fc -matrices 
A^ v ' as PEPS projectors. Acceptable probability of failure e. 
Output: PEPS \ip) with probability at least 1 — e. 

1. t <- 

2. \tpt) entangled pair for each edge e £ E. 

3. H t = E e6B ft 

4. For each v £ V according to total order: 

(a) Ht+x <- H t 

(b) For each neighbor v £ V of v: 

• remove term H^"' v ' from Ht+i 

• compute parent Hamiltonian term H^_'^ ^ 
using A^l 

• add term flfe" to Ht+i 

(c) AddtfW =c(l-pW)toH t+ i 

(d) | Vt+i) measure H t +\ on \ip t ) 

(e) While measured energy nonzero: 

i. \ij>t } 4— measure H t on IV't+i) 

ii. IV't+i) *~ measure H t +i on |"0t ) 

(f ) t 4- t + 1 

FIG. 1. Algorithm constructing injective PEPS 

tonian Hq consisting purely of terms projecting onto 
H e = 1 — ^ j=i f° r each edge of the interac- 

tion graph (step 3). Note, that this simple Hamiltonian 
is gapped. 

We now describe the main iteration of the algorithm 
(step 4), which is illustrated in figure 2. In steps 4a-4c, 
after having selected the next vertex v of the interac- 
tion graph according to the total vertex order, we con- 
struct a new Hamiltonian H t +i from H t : First, we se- 
lect a d-dimensional "physical" subspace from the D k 
dimensional space at each vertex v. This subspace is 
represented by projector P^ y - Then we remove for each 

neighboring vertex v' of v the term Hj: v ' v ' . These are ei- 
ther trivial H e terms or temporary boundary terms (see 
below). Next we compute the new parent Hamiltonian 

terms ^ according to [ ] and add them to H t +\ re- 

flecting the application of A^ . Restricted to the "physi- 
cal" subspace P!$ y , each H^ v,v ) is simply a sum of 2-local 
terms over all edges e from v to vertices v'. Note, that 
parent Hamiltonian terms H^ v ' v ' towards any open "vir- 
tual index" v' are only temporary boundary terms which 
are computed in exactly the same way just as those for 
any other vertex by assuming the identity as the appli- 
cable PEPS map. Since the identity is trivially invert- 
ible, each intermediate PEPS is also injective and thus 
the unique ground state of the intermediate Hamitonian 
H t+ i. Since the "physical" e?-dimensional space is just a 
subspace of the D k "virtual" space that is in fact also im- 
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FIG. 2. In each step, the algorithm processes one vertex v. Ht is grown into Ht+i by removing all existing terms referring 
to v, before the 2fc-local parent Hamiltonian terms are added implementing PEPS projector A^"' at v. Note, that terms around 
v connecting to an open "virtual" index v' (bonds with dotted border) are only temporary and are removed in later steps 
of the algorithm. All terms constraining a vertex v only restrict the physical subspace Pp^ y , while degrees of freedom in the 
orthogonal subspace (1 — ) are eliminated with an additional penalty term that is added per vertex. 



plemented physically in this algorithm, H^ v ' is actually a 
sum of 2fc-local projectors. In order to ensure we produce 
a state with a single d-dimensional local space associated 
to each vertex v in the final PEPS, we add an extra term 
if 



(v) _ 
phy 



— Pphy) in this step. This term penalizes 
the orthogonal complement of the chosen subspace with 
some energy c> A. 

Note, that prior to the execution of step 4d, the sys- 
tem is in the ground state \ipt) of H t by construction. 
This ground state is unique by the injectivity assump- 
tion we make for each intermediate PEPS \4>t) prepared 
in each iteration. In order to transition to the ground 
state \tpt+i) of Ht+i, we run the phase estimation [14] 
algorithm for Hamiltonian H t+ i, perform a binary mea- 
surement to project \ipt) onto the zero/non-zero energy 
subspaces of H t +\, and uncompute the phase estima- 
tion (step 4d). This step requires an inverse eigenvalue 
gap A -1 between these two subspaces that scales with 
0(poly(\V\)) for the phase estimation to be efficient and 
precise enough [15]. We assume that such a gap exists 
for each intermediate parent Hamiltonian H t that we con- 
struct according to the total vertex order defined on the 
interaction graph. 

If the measurement results in the projection onto the 
zero-energy subspace of H t +\ we proceed to the next it- 
eration (step 4e). By Lemma 2, this event occurs with 
probability at least k(A^)~ 2 , where k(A^ v ') is the con- 
dition number of PEPS projector A^ associated with 
vertex v. Note, that the injectivity property of the PEPS 
assures, that each k(A W) is a positive constant. If the 
measurement projects onto the excited subspace of -fft+i, 
we undo the measurement by measuring H t again (step 
4(e)i). If this second measurement results in a projec- 
tion on the ground state, we have exactly undone the 
(unsuccessful) measurement of -Ht+i, otherwise the sys- 
tem is in the excited subspace of H t . In both cases the 
projection onto the groundstate of H t +i can now be at- 



tempted again, with success probabilities k(A^)~ 2 and 
1 — n{A^)~ 2 , respectively (step 4(c)ii). By Lemma 3, 
the inner loop will succeed in projecting onto the ground 
state of H t +\ with probability at least 1 — ^ after at 



Once all 
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most k 2 s attempts, with s chosen as s = j^. 
|V| vertices have been covered, the outer loop terminates 
with the PEPS \tp) in its output register with probability 
at least 1 — e, as shown in Theorem 1. 

Bounding the transition probabilities. — As a first step 
in our analysis, we need a lower bound on the transition 
probability from \ip t ) to \tp t +i)- To this end we proof the 
following lemma. 

Lemma 2. Let \ip t ) — ^TgH^t) ^ e ^ e normalized PEPS 
\A t ), where \A t ) is the unnormalized partial PEPS re- 
sulting from the contraction of PEPS projectors A^ ' for 
all vertices v processed in the algorithm up to step t 



and let Z f 



(A t \A t 



Let \A 



t+ii 



A t +i\A t 



A t+ i is the PEPS projector of time step t + 1. 



l(^+i|^)| 2 > 



> 0. 



Proof. A simple calculation shows 
1 1 



(ipt+i\ipt) = 



<^t|4fil^t 



> 



1 i (MAl 



iA+i\A t ) 



rzt \/Zt+i 

1 1 



Zt+i 



^Zt \/ Z t +i 

1 



c(^+l) 



where 
Then 



(1) 

(2) 
(3) 



(4) 



^max 

where the inequality follows from the operator inequali- 
ties A t+l > and At+ ^ < 1. This implies 



i(^+iiv> t )r > 



i 



J t+X 



t+X) 



(5) 
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But 

Z t+1 = (A t \A 2 t+1 \A t ) > a min (A t+1 ) 2 (A t \A t ) (6) 
= a min (A t+1 ) 2 Z t . (7) 

Thus Eq. 5 and Eq. 7 yield the claim 

p= > ( a ^ t+l) X = (8) 

V^max^t+l)/ K{A t+i y 

Finally, the injectivity assumption of PEPS \4>t+i) im- 
plies left invertibility of A t+ i for each v, thus n(A t+ i) is 
finite, therefore p > 0. □ 



Bounding the convergence rate. — In this section we an- 
alyze the termination probability of the loop at step 4c. 

Lemma 3. Let H t ,H t+ i be Hamiltonians with unique 
zero-energy groundstates \ipt) and respectively. 
Let s be a positive integer. If the system is in state 
\ipt) initially, then the measurement process alternatingly 
measuring H t +i and H t and stopping once \ipt+l) * s 
reached, takes the system to state \ifit+i) with probabil- 
ity at least 1 — after at most s/p alternations, where 

Proof. Let P, Q be the ground state projectors of H t and 
H t+ i, respectively, and let P^ = 1- P, Q 1 - = 1 - Q. 
By Jordan's Lemma, there exists an orthonormal ba- 
sis in which the Hilbert space decomposes into (1) two- 
dimensional subspaces Si invariant under both, P and 
Q, and (2) one-dimensional subspaces Tj on which PQ 
is either an identity- or zero-projector [ ]. 

Since we know that \tp t ) and \ipt+x) are the unique 1- 
eigenstates of P and Q with overlap y/p, exactly one Si is 
relevant to our analysis. This two-dimensional subspace 
is spanned by both, \ip t ) and some |i/'t L )> as well as by 
\ipt+\) and some IV't+i)- Among these four vectors, we 
have got the following relationships [9]: 

Wt) = -VP\M + V^PHt+i) (9) 

iv t x ) = yr^piv^+i) + Vf#m> (io) 

l^+i) = -Vp\i>t) + V^p\^) (ii) 

= + y/ftti) (I 2 ) 

Considering these symmetrical relations, we see that 
alternating measurements of H t and Ht+i generate a 
Markov process among these four states. Since the pro- 
cess terminates whenever it hits j^+i), the only histories 
which can keep the process from terminating are those 
with an initial transition \ip t ) — > and which then 

keep repeating either one of the following two pairs of 



transitions 

-+ -> \4>t +1 ) (13) 

-> |^) -> (14) 

which occur with probabilities (1 — p) 2 and p 2 , respec- 
tively. Thus the process terminates after at most 2m + 1 
measurements with probability 

Ptcrm(p, m) = 1 - (1 - P )(p 2 + (1 - p) 2 ) m . (15) 

To lower-bound this probability we upper- 
bound Pfail(p,™) = 1 - Pterin (p, m) aS 

Pfni\(p,m) < (1 — p)exp{— 2mp(l — p)) which follows 
from (1 - q) m < e-« m , for < q < 1 and m > 0. 
Finally we choose m as a multiple of 1 and find 

Pfaii(p, s/p) < 2^j, which can be seen by straightforward 
calculus. □ 

Proof of Theorem 1. We complete the proof of Theorem 

1 by using Lemma 3 for bounding the failure probability 
Pfaii of the inner loop to derive a lower bound on the suc- 
cess probability of the outer loop over all vertices in V. 
That is, we have to show that (1 — Pf a ii)' y ' > 1 — e. Since 
(1 - Pfaii) |y| > 1 - | Vjjpf a n by truncating higher-order 
terms from the binomial series and assuming |V| > 1 it 
suffices to show |V|pfaii < s. Using Lemma 3, we find the 
first inequality of |V"|pf a ;i < ^ < e, while the second in- 
equality is satisfied by choosing s > Thus, for the al- 
gorithm to succeed with at least probability 1 — e we have 
to choose m > - > p-^-. Since we know from Lemma 

— p — 2pee 

2 that p > -^72, choosing m > > 2pee sumces - 

Thus the inner loop performs at most 2m + 1 < 
measurements. The outer loop iterates over |V| ver- 
tices, thus the total number of measurements is less than 
— ^— + \V\. Bookkeeping of the active Hamiltonian 
terms in the outer loop requires a total time of 0(|V|fc) 
using simple arrays as data structures, and 0(\V\kd 6 ) 
to compute all parent Hamiltonian terms, both of which 
are dominated by the 0(|U| 2 ) time of the inner loop for 
small d. Finally, since each phase estimation step re- 
quires 0(\E\ 2 /A) [10, 14, 15], where O(-) suppresses more 
slowly growing factors such as ea;p(y / ln(|£ , |/A)) [16], we 

find a total runtime of Q( ^§^ + \V\kd 6 ). This com- 
pletes the proof of Theorem 1. □ 

Conclusion. — In this Letter we have shown how to con- 
struct quantum states described by injective PEPS in 
polynomial time by first reducing the problem to the gen- 
eration of a sequence of unique ground states of certain 
Hamiltonians and then preparing that sequence. In fu- 
ture work we will focus on extending the class of prepara- 
blc PEPS and possible performance improvements follow- 
ing from the results of [18-20]. 

This work has been supported by Austrian SFB project 
FoQuS F4014. 
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